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Abstract 

A large number of problems in production, planning and 
scheduling, location, transportation, finance and engineering 
design may be cast into a linear programming framework. 
Inconsistency of constraints and unboundedness of the 
feasible set are of common occurrence while solving these 
problems. This is due to technological environmental and 
competitive factors which interact in a complicated fashion. 
While anti-cycling rules are implemented in most software 
the most used technique for solving linear programs (the 
Simplex method), nothing has been done against the 
vacuousness or the unboundedness of the feasible set. In this 
paper, we discuss a pragmatic approach for progressing 
toward a satisfying solution in line with the decision maker's 
temperament and system of values. At the heart of this 
approach lies in the use of the fuzzy paradigm to represent 
Decision maker's views and to handle them. Numerical 
examples are carried out for the sake of illustration. 
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Introduction 

Linear programming is an important tool in the 
arsenal of means at a Decision Marker's disposal. 
Indeed it is one of the most frequently applied 
Operational research modals in solving real word 
problems [(R. J Vanderbei, 2008), (S. Nash, 1996), (G. 
Dantzig, 1951)]. Linear programming's theoretical 
underpinning is now well established and as a result, 
a broader array of techniques including the Simplex 
method (S. I Gass, 2003), the Ellipsoid method (D. 
Dubois, 2011) and the Interior point methods (N. 
Karmarkar, 1984) have been developed. User-friendly 
software with powerful computational and 
visualisation capabilities has also been pushed 
forward. While solving a real-life problem using 
linear programming, it may happen that either the 
constraints of the problem are inconsistent or the 
feasible set is not bounded. These irregularities pose 
serious challenges, to the Decision Maker, mostly 
when he is in need of a solution "here and now". 
While anti-cycling rules serve as guidelines in 
implementing efficient techniques for dealing with 



degeneracy, there is notable lack of operational tools 
for guiding decision making in the presence of 
inconsistencies and unboundedness. 

In this paper we take the option of facing these 
irregularities by abandoning the pure rationality 
principle (A. Schuts, 1943) and its naturalist trend of 
seeking an "optimum optimorum", that anyway does 
not exist in these cases. 

Rather, we adopt a constructivist perspective of 
creating a satisfying solution (H. A Simon, 1997) in 
line with preferences and system value of the 
Decision Maker. At the heart of our approach lies in 
the fuzzy paradigm (L. A Zadeh, 1978) that helps in 
representing and handling Decision Maker's 
preferences. 

For other uses of Fuzzy set Theory in Decision Science 
we refer the reader to [(L. G Khachiyan, 1979), (S. A 
Orlovsky, 1978)]. 

The remainder of this paper unfolds as follows. In the 
next section we briefly discuss how to check 
infeasibility and unboundedness for a linear program. 
In section 3 we present some ways for guiding 
decision making in the presence of inconsistency and 
unboundedness in a linear programming setting. 
Numerical examples aiming to illustrate the 
effectiveness of our approach are carried out in 
Section 4. We end up in Section 5 with some 
concluding remarks along with lines for further 
development, in the field. 

Checking Infeasibility and Unboundedness in 
a Linear Programming Context 

Before proceeding further, we find it convenient to 
discuss some ways for checking inconsistency of the 
constraints and unboundedness of the feasible set, of a 
linear program. An interested reader may consult [8] 
for more details on these matters. 

Inconsistency Criterion 

The following result due to Minkowski and Farkas (S. 
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Kutateladze, 2010) and bearing the latter's name will 
be of paramount importance in stating Inconsistency 
criterion that will be used in this paper. 

1) Farkas Lemma 

Let A be a given nxm real matrix and b 
given m vector. The inequality b T y > holds 
for all }g1° satisfying Ay > if and only if 
there is an n vector 8 > such that A T 8 = b. 

The proof of this result may be found elsewhere 
(M. Comforti, 2008). 

A result that extends quite canonically the Farkas 
Lemma is known under the name of S-Lemma (I. 
Polik, 2007). The basic idea of this widely used 
mathematical result came from Control Theory but 
it has important consequences in quadratic and 
Semi-definite Optimization. 

It is worthy nothing that Farka's Lemma is 
equivalent to the following statement: 

The system 

[ A T x = h 

(1.1) 

is incompatible if and only if the system 
b T y<0 




Ay>0 

y e R m 



(1.2) 



has a solution. 



2) Feasibility Test for a Linear Program 
Consider the following linear program 



(P) 



min cx 
Ax = b 
x>0 



It is not a loss of generality to ponder only equality 
constraints. As a matter of fact, by making use of 
slack variables, a system of constraint inequalities 
can be converted to a system of equalities. 

From Section 2.1.1 we can derive the following 
inconsistency criterion for (P). 

Solve the system (1.2). 

If this system has a solution then (P) is 



inconsistent. Otherwise (P) is feasible. 
Unboundedness Criterion 

The following result due to Weyl tells us that the 
unboundedness case may occur. 

1) Weyl's Theorem 

Consider X = [x e M" / Ax < b,x > 0} 

where A and b 

are as in Section 2.1.1. Then either X is empty or 
X is a bounded polyedra or X is an 
unbounded polyedra set. 

For the proof of this result the reader may consult 
[12] 

2) Case of Non-finite Optimal Solution 

Consider (P) and let B be an mxm 
submatrix of A such that det 5^0 i.e. B is a 
basis. Let 1(B) stands for the set of indexes of 
columns of B and 



iel(B) 



Denote by k the index of the basic variable that 
enters the basis i.e k is such that 



Z k ~c k = max (Zj-Cj) 



(1.3) 



Then the following statements, the proof of which 
may be found in (S. Nash, 1996) hold true. 

Proposition 1 

If Zj —Cj<0 Vj, then the basic solution 
corresponding to the basis B is optimal. 
Proposition 2 

If for k defined in (1.3), a tt < Vi then (P) 
has no finite optimal solution. 

Dealing With Inconsistency and 
Unboundedness in a Linear Programming 
Framework 

In an attempt to apply Linear Programming to 
concrete real-life problems, one may (according to the 
Weyl's Theorem) encounter situations in which the 
feasible set of the linear program under consideration 
is either empty or not boundless. In this case the 
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Decision maker will not be happy with inconsistency 
and unboundedness reports. Our claim here is that 
simple and intuitive ideas from Fuzzy set Theory may 
be used, with proper reason; to make a relaxation of 
constraints in the case the feasible set is empty. The 
same ideas may also help convert the original 
problem into a flexible one, in the case of 
unboundedness. Both approaches rely heavily on 
Decision Maker's views and opinions. In this way one 
enhances the degree of conformity between the 
evolutions of the decision-making process and the 
value system of those involved in the process. The 
ultimate objective is to deliver a satisfying solution 
instead of an optimal one, which anyway does not 
exist. 

Resolving Inconsistencies among Constraints 
1) Analysis 

Consider the linear program 



(P2) 



min cx 

Ax<b 
x>0 



where A is an nxm real matrix, b an n 

vector, c an m — vector and x e R m . Suppose 
that the feasibility test (see Section 2.1.2 ) reveals 
that the feasible set of (PI) is empty. Instead of 
sending an inconsistency report, one may attempt 
relaxations of the objective and the constraints in a 
way to find a satisfying solution. Such a relaxation 
is given as the following Flexible program: 



(PI)' 



min cx 

Ax^b 
x>Q 



Where ~ means that min and < are no longer 
strict imperatives but some leeways may be 
accepted. An interpretation which goes with 
(PI)' is as follows. 

Find XEl-jl6R m ,l>0j such that cx is 
as 

well as possible below a reasonable level b {) and 
such that the constraints 



A ; x<^ ( .; i = \,...,n 



are met as well as possible. Or merely, find x & X 
such that the following system holds 



AjX^br, r = 0,l,...,n 
Where Aq = c. 



(1,4)' 



l _A i x-b i 

d, 



A convenient way to represent these soft 
constraints is through appropriate fuzzy sets of 
M, the membership functions of which are 

i = 0, 1, n. 

These membership functions may be defined as 
follows: 

if A,.x < h 

(i.e if the ("constraint of (1.4) is satisfied), 
if b i < A.x<h +d { 

ju, (x) = -J (i.e if the /"'constraint of (1 .4) 

is violated to some acceptable extent) , 
if AfX > b. t + d ; 
(i.e if the /"'constraint of (1.4) is strongly violated). 



Where d i z = 0,l,...,n are fixed by the Decision 

Maker. They are subjectively chosen constants of 
admissible violations. 

According to the Bellman-Zadeh's confluence 
principle (R. E Bellman L), a decision in a fuzzy 
environment is an option that has the highest 
membership degree in the fuzzy set intersection of 
the objective and the constraints of the problem. 

This means we seek for x* = arg max jU D (x), 



Where jU D (x) = A jU^x), A standing for the 

i=0 

connective 

"and". Using the min 7 -norm (L. A Zadeh, 
1978) to translate the meaning of the intersection, 
we have to solve the following optimization 
problem 



(1.4) 



(P2) 



maxmin ;=0 n //,(x) 

n 

x e f^SuppjUi 



i=0 
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Where Supp/U i stands for support of jU r 
(P2) 

v ' is equivalent to the following mathematical 
program: 

max A 



(P2)' 



2<1_A^A;/ =0 ,1,..., W 



Ax<b t +d i ;i = 0,l,-,n 
x>0 



This linear program can be solved by existing 
mathematical software like LINDO, LINGO, 
COMPLEX, XPRESS. 

2) Description of the Method 

From the foregoing discussion we obtain the 
following flowchart for dealing with inconsistency 
in a linear programming setting. 



Start 




Read A, B, C 



Solve the system (1.2) 




Solve the LP using LINGO, 
CPLEX or XPRESS 



1 


T 


Print the solution 







(^Stop^) 




Put (PI) in the form of (PI)' 



Ask the DM to fix 

d ;z = 0,1,. ..,m 



Ask DM to fix 

d' > d t ;i = 0,...,n 





Print the solution 




r 




Stop 



FIGURE 1 FLOWCHART OF LP METHOD IMPLEMENTING 
INFEASIBILITY 

Handling Unboundedness of the Feasible Set 
Consider the following linear program 
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(P3) 



max cx 

Ax<b 

x>0 



functions. 



where A is an n x m real matrix, b an n — vector, 

c an m— vector and x e IR m . Suppose that the 
criterion for detecting unboundedness ( see Section 
2.2.2) reveals that (P3) has no finite optimal 
solution. Instead of delivering an unboundedness 
report, one may, with the help of the Decision Maker, 
reformulate (P3) as follows: 



Find xeX={xeM. m /x>0} 
such that cx is sufficiently large 

while Ax < b is satisfied as well as possible 



ju (x) 



cx-(b Q -d ) 



if cx < b - d 
(i.e if cx is too small), 



if b - d <cx<b 



1 



(i.e if cx is enough close to b ), 
if cx > b n 



(i.e if cx is sufficiently large). 



and for i = l,...,n : 
1 



J u i (x) = 



if Ax < h 



\- AiX bi if b i <A i x<b i 



if A ( jc > b t + d i 



Or merely as the following Flexible program: 



(P3)' 



max cx 

Ax^b 

x>Q 



Following Zimmerman (H. J Zimmermann, 1996), 
(P3)' may be written as 



(P3)" 



Find xe X such that 

cx ^ b Q 

A t x ^ b{, i = l,...,n 
x>0 



Where b is a threshold provided by the Decision 

maker to calibrate his view on the fact that the 
objective function is large enough. These flexible 
constraints may be represented as fuzzy sets of X 
with membership functions: 

These membership functions should be determined by 
capturing the Decision maker's preferences. 

For the sake of illustration assume that the decision 
maker agrees with the following piecewise linear 



Start 




ReadA,B,C 



Determine a basis B 
along with 1(B) 



1 

Find the corre 
Solution x B 


sponding basic 


► 


T 






FIGURE 2 PROCEDURE FOR IMPLEMENTING 
UNBOUNDEDNESS IN A LP PROBLEM 
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1 




Find the index k of the 
entering variable 







Find the "basic solution at hand 
Is optimal" 





Stop 



Find the index r of the 
outgoing 




Make simplex transformation 
to obtain a new basis B' 



Print solution 
of the problem 
is infinite 



Find the corresponding basic 
solution x R , 



Put B = B' and go to s — \ 



FIGURE 3 PROCEDURE FOR IMPLEMENTING 
UNBOUNDEDNESS IN A LP PROBLEM CONTINUED 



x = arg max KeX ju D (x) where 



jU D (x)= min 

i'=0, !,...,« 



Section 3.1.1 and translating the intersection by the 
min operator, the solution to (P3) will be 

Such an x may be obtained by solving the following 
mathematical program: 



(P3) m 



max min, =0 !„//,(*) 

n 

xef^Suppu 



i=0 



(P3)'" is equivalent to the following optimization 
problem 

max A 

x < cx-(b -d ) 
d 

(P3f l- A ' X ~ bi ; i = l,...,n 

d i 

cx> b -d o 
AjX < b i +d t ; i = l,...,n 
^;>0 

The above considerations lead to the following 
procedure for implementing unboundedness in a 
linear programming framework. 

Numerical Examples 

Case of Inconsistent Constraints 
Consider the linear program 

max (5xj +4x 2 ) 



OP/) 



Xj + x 2 < 2 
-2x i - 2x 2 < -9 
Xj > 0, x 2 > 



This program is infeasible. Indeed, the second 
constraints implies that 



Applying Bellman-Zadeh's confluence principle like in 



x x + x 2 > 4.5 



Which contradicts the first constraints. Let now apply 
the method described in section 3.1 to reconcile theses 
inconsistent constraints. Assume that in the context 
under consideration the Decision Maker interprets 

max (5x, +4x 2 ) as 5x, +4x 2 should be, as well as 
possible, greater than 10. 
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Then reads: 

Find ieX=j(i p i 2 )6K 2 /i | >0,l 2 >0} such 
that: 

5x, + 4x 2 > 10 (0) 

Xj + x 2 < 2 (1) 

-2xj-2x 2 <-9 (2) 

Suppose further that the Decision Maker fix the 
respective constants for admissible violations of 
constraints (0) — (2) 

As follows d =5, d x =\, d 2 =4.5. 

Then the flexible constraints (0) — (2) may be 
represented as fuzzy sets of X the membership 
functions of which are respectively: 



JU Q (x) = 



1 if 5jc, + 4jc 2 > 10 

5x,+4x 2 -5 if5 < 5 ^ +4 ^ 2<10 



10-5 



ju l (x) = 



1 



if 5xj + 4x 2 < 5 



if x x + x 2 < 2 



(x,+x~) — 2 . _ „ 
1— ^ if2<x t + x 2 <3 



(2 + l)-2 



if Xj + x 2 > 3 



1 

1- 



-(2^ + 2x 2 ) + 9 



if -2x, -2x 2 <-9 

if -9<-2x, -2x 2 <-4.5 



-9 + 4.5 + 9 
if -2x; -2x 2 >-4.5 



It is clear that the supports of the three fuzzy sets are 
respectively: 



Supp (ju ) = [x e X 1 5xj + 4x 2 > 5} 

Supp (ju 1 ) = {xe X / x { +x 2 <3) 

Supp (ju 2 ) = {xeX / -(2x, + 2x 2 ) < -4.5} 

This resulting program (see Section 3.1) is: 



max X 

^ < 5x t + 4x 2 - 5 
5 

X< l-(x 1 +x 2 -2) 

(-2xj - 2x 2 + 9 



4.5 



X< 1 — 

5xj + 4x 2 > 5 
Xj + x 2 < 3 
Xj + x 2 > 2.25 
x x >0,x 2 >0,2>0 

This program can be merely written: 
max X 

^-X!-0.8x 2 <-l 
X + Xj + x 2 < 3 

l-0.44x!-0.44x 2 <-l 



- 5x, - 4x n 



Xj + x 2 



<-5 
<3 

< -2.25 



vj yv 2 

2 > 0, x x > 0, x 2 > 
Which yields using LINDO software the solution: 



x,= 2.777778 
x 2 = 0.000000 
X = 0.222222 



As can be seen, our approach has reconciled 
inconsistent constraints to deliver a satisfying solution 

X* = (2.777778, 0.000000) and the maximum 
degree of overall satisfaction of constraints is 22%. 

This degree may be improved by a less stringent 
choice of 

d i (i = 0, 1, 2) by the Decision Maker. 

Case of Unboundedness of the Feasible Set 
Consider now the following linear program 
max (Xj -4x 2 ) 



OP) 



-2xj +x 2 < -1 
-x y - 2x 2 < -2 
^ > 0, x 2 > 
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In this program we can set x 2 to zero and x, 

arbitrary large. As long as X, is greater than 2, the 

solution will be feasibly and it gets large, as the 
objective function does too. Hence, the solution of this 
linear program is unbounded. We are going to apply 
the method described in section 3.2 in a way to find a 
satisfying solution in line with Decision Maker's 
preferences. The relaxation program associated with 



(P u ) ^ 



max (x l - 4x 2 ) 
-2xj +x 2 ^ -1 



-Xj - 2x 2 ^ -2 
x l > 0, x 2 > 



or 



(PJ 



Find x e X = {x = (x 1 ,x 2 )/x 1 > 0, x 2 > 0} 
such that 

(x 1 -4x 2 )>& (0) 
-2x!+x 2 <-l (1) 
-x, - 2x 2 < -2 (2) 



Assume ^ () is fixed to 20 by the Decision Maker. 

Suppose also that subjectively fixed constants for 
admissible violation 

corresponding to flexible constraints (0), (1), and (2) 
are 

d = 10, d x = 2 and d 2 = 4 respectively. Then the 

membership functions associated to (0), (1) and (2) are 
respectively: 



jU (x) = 



1 if (x!-4x 2 )>20 

(x!-4x 2 )-10 . 



20-10 







if 10<X!-4x 2 <20 



if (x 1 -4x 2 )<10 



if -2xj +x 2 < -1 



, 1 _^x L +x 2 ) + l if-1< _ 2 <j 
(-l + 2) + l 1 2 







if (-2x 1 + x 2 )>l 



ju 2 (x) = 



1 if - Xj - 2x 2 < -2 

(-x-2x 2 ) + 2 



1- 



(-2 + 4) + 2 



if-2<-X!-2x 2 <2 



if (-\x l -2x 2 )>2 

The counterpart of (P3)" for this example is: 
max A 

A - 0.1^+0 .4x 2 <-l 



X-x x +0.5x 2 



<0.5 



A-<Zi.25x x -<Z).5x 2 <§5 
x l -4x 2 >10 
<-l 
< -2 

A <1 
X > 0, x, > 0, x 2 > 



Using LINDO software, we obtain the satisfying 
solution x* = ( 20, ) with the maximum degree of 
overall satisfaction of constraints of 100%, i.e. A* =1. 

Concluding Remarks 

Methods for solving linear programming problems 
should not focus only on regularities. They should also 
be able to provide guidance on those features that 
contradict regularities and create complexity, like 
inconsistency of constraints and unboundedness of the 
feasible set. Theses singular features pose challenging 
questions to the Decision Maker, mostly when he is in 
need of implementing an action here and now. Dealing 
with inconsistency and unboundedness in a linear 
programming setting may be regarded as a leap in the 
dark. Nevertheless there are different ways to leap in 
the dark. In this paper, we have adopted a 
constructivist position to obtain a satisfying solution in 
a turbulent environment where either inconsistency or 
unboundedness is in the state of affairs. The grim 
preamble to the construction of such a satisfying 
solution is the capture of the Decision Makers views 
and opinions. Fuzzy sets theory has been used both as 
a language for representing Decision Maker's 
subjective opinions and as a tool for combining them. 
Extensions of this work can explore various directions. 
It is worth mentioning that using the min operator to 
translate the meaning of the intersection may be a too 
pessimistic approach. In decision problems where such 
an ultra pessimist attitude cannot be accepted, 
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compensatory operators can be used. For instance jU D 
can be defined as follows. 



ju D (x) = y min ju i (x) + (I- /)rmn 



Where y £ (0,1) is a coefficient of compensation [14]. 
Moreover, other membership functions that suit better 
the situation at hand may be considered in lieu and 
place of piecewise linear functions described here. The 
most used membership functions include hyperbolic, 
hyperbolic inverse, logistics and cubic functions (H. 
Rommelfangger, 1996). A related problem, important 
but much more difficult is to ponder distances to 
inconsistency and unboundedness (M. J Canovas, 2006) 
in a way to advise the Decision maker's during the 
decision process. This problem deserves more 
thoughts and further consideration. This will become a 
topic for future research. 
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